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Abstract 

We present in this paper a proof of well-posedness and convergence 
for the parahel Schwarz Waveform Relaxation Algorithm adapted to 
an N-dimensional semilinear heat equation. Since the equation we 
study is an evolution one, each subproblem at each step has its own 
local existence time, we then determine a common existence time for 
every problem in any subdomain at any step. We also introduce a 
new technique: Exponential Decay Error Estimates, to prove the con- 
vergence of the Schwarz Methods, with multisubdomains, and then 
apply it to our problem. 

1 Introduction 

In the pioneer work [16], [17], [18], P. L. Lions laid the foundations of the 
modern theory of Schwarz algorithms. He also proposed to use the Schwarz 
alternating method for evolution equations, and studied the algorithm for 
nonlinear monotone problems. Later, Schwarz waveform relaxation algo- 
rithms, by refering to the paper [2j, were designed independently in [12j and 
for the linear advection-diffusion equation. They try to solve, on a given 
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time interval, a sequence of Cauchy Problems with the transmission condi- 
tions of Cauchy type on overlapping subdomains. The algorithm is well-posed 
with some compatibility conditions. 

An extension to the nonlinear reaction-diffusion equation in dimension 1 
was considered in [10]. For nonlinear problems, especially evolutional equa- 
tions, there are some cases that the solutions are blowed up in finite time, 
which means that if we devide the domain into several subdomains, at each 
step we can get different existence times in different domains, and we do not 
know if there exists a common existence time for all iterations. However, 
with the hypothese /'(c) < C in [10], we do not encounter this difficulty and 
the iterations are defined naturally on an unbounded time interval. Proofs of 
linear convergence on unbounded time domains, and superlinear convergence 
on finite time intervals were then given in case of n subdomains, based on 
some explicit computations on the linearized equations. Another extension 
to monotone nonlinear PDEs in higher dimensions were considered by Lui in 
[19] . [20] . In these papers, some monotone iterations for Schwarz methods are 
defined in order to get the convergence of the algorithm, based on the idea of 
the sub-super solutions method in PDEs and no explosion is considered. Re- 
cently, an extension to Systems of Semilinear Reaction-Diffusion Equations 
was investigated in [5]. Some systems in dimension 1 were considered and 
the proofs of the well-posedness and the convergence of the algorithm used 
in this paper were a development of the technique used in [TU] . 

We consider here the semilinear heat equation fl2.ip . in a domain Q = 
D X (a, b) of M^, with the nonlinearity of the form (12. 2p . which allows explo- 
sion of solutions in finite time. We cut the domain into bands Qi = Dx[ai, hi), 
with ai = a and bi = b. These bands are ovelapping, i.e. for all i G {1, / — I}, 
ttj+i < bi < In each of these subdomains, we solve a heat equation with 
Dirichlet limit conditions. Since the domains are not smooth, we cannot 
use classical results about semilinear heat equations on smooth domains; we 
then establish some new proofs of existence for a general domain in Theo- 
rem 2.1. Applying the results in Theorem 2.1 for the equation (12. 7p . we get 
an existence result for a semilinear heat equation in a domain of the type 
n = Dx {a,b) in Theorem 2.2. 

Theorem 2.3 confirms that the algorithm is well-posed and there exists 
a common existence time for all subdomains at all iterations despite of the 
phenomenon of explosion. The common existence time T* is computed ex- 
plicitely so that one can use it in numerical simulations. This is a coUolary 
of Theorem 2.2. 
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We prove in Theorem 2.4 that the algorithm converges hnearly. There 
are five main techniques to prove the convergence of Domain Decomposition 
Algorithms and they are: Orthorgonal Projection used for a linear Laplace 
equation (see tl6j), Fourier and Laplace Transforms used for linear equations 
(see [8], [9], [11], [E]), Maximum Priciple used for linear equations (see [T3]). 
Energy Estimates used for nonoverlapping algorithms (see pQ), and Mono- 
tone Iterations (see [T^], [20]) used for nonlinear monotone problems. The 
convergence problem of overlapping algorithms for nonlinear equations is still 
open up to now. In this paper, we introduce a new technique: Exponential 
Decay Error Estimates, basing on the idea of constructing some Controlling 
Functions, that allows us to prove that the algorithm converges linearly. 



2 The main results 

We consider the semilinear heat equation 

dtu- Au- f{u) = 0, (2.1) 

with the assumptions on /: 
/ is in C^{R) and there exists C/ > 0, p > 1 such that |/'(a;)| < 

(2.2) 

We first set an existence theorem for the initial boundary value problem, 
and more important, new estimates on the solution. We need here some 
notations. We set pi = ^^^^^^ , a = ^i:^ ~" ^)- ^2 are positive 

numbers such that + = ^ aiid hPi < 1- We denote by ll^Hfc,/! the norm 
\W\\l''{o,t,l'^{u)), when u belongs to L^{0,T, L^{u)), where u is some domain 
in M^, k, h can be 00. We define 

r(r,m) = [(47r)-P^ C/ max(l, 2P-2)(4r + 4m) ]-^^ , 

G{r; T, m„m,) = ( ^^^^'"It""'' ) Jo [^f max{l, 2^-'}{m, + c'^)C^ + m,] ''^ dC 

(2.3) 

Consider the problem 

dtw - Aw = f{w + v) inOx(0,r), 
w = ondOx[0,T], (2.4) 

w{.,0) = in 0. 
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Theorem 2.1. Let O be a bounded domain in M^, m{0) is its measure. 
Suppose V G C([0, To, L^](fi)) and \v\ < M a.e.. Denote Ri = 2max|^|<A^ 

P+3 

\f{Q\m{0)2 -g^ny. Then, there exists a local time T^,, = min(To, 1, t{Ri, m{0) ))^ 

4p 

such that for all T < T^, equation ^2.4\) has a unique solution w in L°°{0 x 
{0,T))nC{[0,T],L^{O))nL\0,T,H^{O)) anddtw e L\0,T,L\O)). Fur- 
thermore, 111^1100,00 < M^,, where 

:= (^)^ [ Cf max(l,2P-2) (G'-i(T,)'^ + MP-'m{Of-^) G-\tS~^ 

+m(C)5maxi^l<M 1/(01 ] 
G-\n) = G-\T,; T, MP-^m{0f-^,m{0)-2 max|^|<M |/(C)I)- 

(2.5) 

p+3 

Suppose V e C([0, To], L2(^]))^L-(0, To, L^^in)) a.e.. Denote R2 = 2^||/(t;)||oo,2. 

Then, there exists a local time T* = min(To, 1, T(i?2, \ \'^\fL°^{o t L2p(n)) > ■^'"^^ 
that for all T < T^: , equation i \2.4^ has a unique solution w in L°° {Ox {0,T))n 
C([0,T],L2(0)) nL2(0,T,ifi(C')) and dtw e L2(0, T, L2(0)). Furthermore, 
||w^||oo,oo < M^, where 

M. := (^)^ [ Cf max(l,2^-2) (^-^(T.)^ + llt;^^-^) ^-^(T.)^ 



+ 11/(^^)1100,2] 
G-HT.) = G-i(T,; T, ||/(^;)|U,2) 



(2.6) 



We consider now a bounded domain of the following form Q = D x 
(a, b) C M^, where D is a bounded domain with smooth enough boundary 
dD in M^^^. The boundary dQ of Q is made of three parts, = D x {a}, 
Tfi = D X {b}, and Tc = dD x [a,b). Dirichlet data g are given on the 
boundary dQ x (0,T), defined by ql on Tl, qr on Vr, gc on Tc- These 
functions are all continuous. We now introduce the basic initial boundary 
value problem for ( 12. 4p : 

dtu -Au = f{u) in X (0, T), 

u = g in<9fix(0,T), (2.7) 

u{., 0) = Uq in fl. 



4 



Theorem 2.2. Let Uq in C{Q) and g in C{dQ) with Uq^qq = g\t=o- Let M be 
a positive constant such that M > max{\\uo\\oQ,\\g\\oo) ■ Let Ri be like in The- 
orem \2.I[ Then, there exists a limit time = min(l, t{R, Mp~ m{il) )), 
such that for all T < T^,, equation l \2. 7p has a solution u in L°°{Q x (0, T)) Pi 
C([0,^],L2(^])) n L\0,T,H ^{n)) and dtu G L\0,T,L\n)). Furthermore u 
is continuous on x (0,T) and \ \u\\oo,oo < M + M* where M* is obtained 
from ( \2. 5p by replacing O by Q. 

We divide the domain Q into / subdomains with Qi = D x (oj, bi), with 
ai = a and 6/ = b. We suppose for each ? G {1, ...,/— 1}, Oj+i < b^ < 
we denote by Lj the length of Qf Li = hi — ai, and by Si the size of the 
overlap Si = hi — Oj+i. 

The parallel Schwarz Waveform Relaxation Algorithm solves / equations in 
/ subdomains instead of solving directly the main problem (12 .7^ . The iterate 
in the j-th domain, denoted by u'j, is defined by 

' dtu'; - Au'; = f{u]) in fi,- X (0, T), 

< u]{-.aj,-) = u)z\{-,a,,-) inDx(0,r), (3.8) 
u]{-,b,,-) = u)-l{-,b,,-) inDx(0,r). 

Each iterate inherits the boundary conditions and the initial values of u: 

u] = g on dQj n dQ x (0, T), u']{-, -, 0) = uq in Qj, 

which imposes a special treatment for the extreme subdomains, 

u'l{-,a,-) = g, u'}{-,b,-) = g. 

An initial guess is provided, i.e. we solve at step equations ( 12. Sp . with 
boundary data on left and right 

M;(-,ajv) = g^ in D X (0, T), 6^-, ■) = /i° on D x (0,r). 

Definition 2.1. The parallel Schwarz waveform relaxation algorithm is well- 
posed if there exists a local time T* < such that for all T < T* , each 
subproblem ( \2. 8h in each iteration has a solution over the time interval (0, T), 
and the set of solutions {u^, j G 1,/, G A^} is bounded in C{Q x [0,T]). 



5 



Let M a positive number. According to theorem \'2.'2\ the following prob- 
lem has a solution (pM in some interval CiVt x [0, Tq]): 

dt(t)M - ^(i)M = f{(t)M) in X (0, To), 
(t)M = M inafix[0,To], (2.9) 

(^m(-,-,0) = M mn. 

The next theorem gives a common existence time for the iterates: 
Theorem 2.3. Let M he a positive constant such that 

M > max(||Mo||oo, IblU, (||5'j°||oo, | 1 U, ^ e 1, /) ). 

Suppose the data uq and gj, gj and are continuous and satisfy the com- 
patibility conditions 

Uo\dQ = 9\t=o, Mo(-,«i) = ^^Jii^^g, "o(-,f'i) = ^?l^^o- 

Let be greater than the maximum of 4>m on the boundaries of Q x (0,T) 
and Qj x (0,T), j e TJ. PutT* = min(To, 1, T^,t{Ri, MP-^m{nj)^),j e 
1, /)). Then the parallel Schwarz waveform relaxation algorithm i\2.8^ is well- 
posed with a local time at least equal to T* . 

We denote by Cj the error Uj — u, where u is the solution of (12. 7p . Let 
P be a function from M to M such that (i) P G C^{R); (i) P{x),P"{x) > 
Vx e M; (iii) P{x) = iff a; = and P'(0) = 0; (iv) VM > 0, there exists 

K{M) such that < K{M), V x G M, |a;| < M. We finally state the 

convergence of the algorithm: 

Theorem 2.4. With the same assumptions as in Theorem \2.S\ let'-/ be a con- 
stant large enough and denote by e the constant vi r ^ /'^ ■ If put Ek = 



2 L2...L1. 

T^^^jGU W^i^j) 6xp(-70llc ffi,x(o,T) ) = "^ax^-gT7 ll^(«j-«)exp(-7^)llc(fv^) 
and Ek = max^- ^ ^ j_^ {Ek+j}, the parallel Schwarz waveform relaxation algo- 
rithm ( \2.8^ converges linearly in the following sense: for every T < T* , 

En < -Eoexp(-ne), Vn G N, 

and as a consequence 
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3 Proof of The Existence Results for Semilin- 
ear Heat Equation in an Arbitrary Bounded 
Domain - Theorem 12.11 



3.1 Premiliary Results 

Let u be an arbitrary bounded domain in and consider the following 
equation 

atC-AC = ina;x(0,T), 

C(.,.) = onaa;x[0,r], (3.1) 
C(.,0) = Co oncD. 

We consider the operator = -^C, D{A) = {C|C e H^{u),AC e 
L'^{u)}. According to Proposition 2.6.1 page 26 [3J, A is an m-dissipative 
operator with dense domain in L'^{u). Let S{t) be the Dirichlet semigroup 
associated with A on L'^{uj) (see [3]). If Co ^ D{A), due to Theorem 3.1.1 
page 33 [3], C{t) = 5'(t)Co is a solution of (13.11) . 

According to Section 6.5 page 334 |6j, there exist a sequence of eigen- 
values < Ai < A2 . . . and Afc — i- 00 as — )• 00 and an othornormal basis 
{wfcjfcli of L'^{uj) which is also an orthogornal basis of Hq{u), where Wk is 
the eigenfunction corressponding to A^: 

Awk = XkWk in u, 
Wk = on du. 

Denote by <, > the scalar product in L'^{u), we have the following Lemma 

Lemma 3.1. If (0 G L'^{uj), then we have that 

00 

S{t)Co = J2 < ^0, Wk > Wk. (3.2) 

k=l 

Proof. We prove the Lemma in two steps. 
Step 1: Suppose that Co ^ D{A) and put 

00 

u{x, t) = e"*^^- < Co, Wk > Wk, 

k=l 
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we prove that is a solution of f l3.ip . 

We compute the norm of jy{-,t) in L?{<^) 

oo oo 
fc=l k=l 

The inequahty imphes G L°°(0,T, L^(a;)). 

We now need that z/ G lv^(0, T, ifQ(a;)), and it is enough to prove 

oo „T 

/ (e"*^' < Co,«'fe >f\\wk\\]ii(^^)dt < oo. 

A:=l 

We estimate the term on the left hand side of the previous inequality 

oo „x oo „x 

k=l "^0 k=l "^0 

oo 

k=l 



The estimate leads to G -L^(0, T, Hq{uj)). 

Since i^(0, .) = YlT=i ^ Co)""^*: > Wk = Coi we only have to prove that 
dtiy — Az/ = 0. In order to do this, we put z/„ = J2k=i(^~^^'' ^ Coi'W^fc >)wk- 
We compute by a direct manner 

n n 

5tM„ = ^-Afc(e"*^* < Co, Wk >)wk = ^(e"*^'= < Co,Wk >)Awk = Az/„. 

k=l k=l 

It follows from the previous equation that 

ni^dtVnfdxdt = — / Vi'n'V{dtiyn)dxdt 
Jo Juj 



The inequality /^(Vz/„)2(0)dx = L ELi < Cc^^fc (Vw;,)^^^ <£ ^^=1 < 
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Co.Wk >^ [Vwkfdx = UVCo^dx = IIColl^w^) implies JJdtUnYdxdt < 



llColl^i(^). We obtain a*// G L\0,T, L\u)) &nd then dti^ = Ai^ in L\0,T, L\u)). 
Step 2: We prove (lO) . 

We now put jji = u — Q and recall that ^ = S{t)C,Q is a solution of (13. ip . It 
follows that /i is the solution of the following equation 

' at/i-A// = inwx (0,T), 

= on aw X [o,r], 

0) = on w, 

[ ^^ e L\0, T, H^,{u)), dtiM G L2(0, T, L2(a;)). 
We derive from fl3.3p that 



(3.3) 



< dtfi{., s), fi{., s) > ds + / < V/i(., s), V/u(., s) > cfs = 0. 

'0 

This means iMtWm^) + lo W"^ f^Whi.) = 0- Then = 0, or u = (; and O 
holds for Co e D(A). Since D(A) is dense in L^{uj) and S{t) G ^(^^(cu)), 
(13:21) holds for Co e L2(a;). □ 

Lemma 3.2. Suppose that g G C([0, T], L^((X')) and po ^ D{A), we consider 
the following equation 

dtp — Ap = g inux(0,T), 

p(.,.) = ondujx[0,T], (3.4) 

p(.,0) = po on u. 



Then /ias a solution m L'^{0,T, H^{u)) n L°°(0, T, /.^(w)) given by 

the following formula 



p{.,t) = S{t)po+ S{t - s)g{.,s)ds. 
Jo 



(3.5) 



Proof. According to Lemma [3?H S{t)pQ is the solution of (13. ip with the initial 
condition po and S{t)po G ^^(0, T, H^{u))nL°^{0, T, L^iu)). We can suppose 
with out loss of generality that po = in our proof. Then, by Lemma I3.lt 
formula (13.51) can be rewritten in the following form 



p{.,t) = Y^ / e-^'"'^'^ <g{.,s),w,>ds 
k=i L-^o 



Wk. 



(3.6) 
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step 1: We prove that p G 1^(0, T, HI{uj)) n L'^lO, T, L^co)). 
The fact that p e L°^(0, T, can be derived from the following 

inequality 



o<^ r 

V / e-^'-^^^^ <g{.,s),Wk>ds 
k=i L-^o 



2 <x, 
< 



^ ^llfl'llz,2(o,r,L2(w))- 



Moreover, from 



°° r ft 



E 



-(t-s)Afc 



< 9{-,s),Wk > ds 







k=l 

oo r /-t 



= I] / e-^*-^)^'^ <^(.,s),«;fe>ds 
fe=i L-^o 
°o /■< /■* 
< 5^ A,. / e-'^'-'^^'ds <g,Wk>^ds 
k=i -^0 

we obtain that p G L^(0, T, 7^0(0;)). 

Step 2: Wc prove that dfP — Ap = g. 

Put pn{t) = Efc=i[/o < 9{;S),Wk > ds]wk and = XlLi < 

g{t),Wk > Wk. Since g G C([0, T], ^^(a;)), 9tp„(i) = Ap„(t) + This 
formula implies 

/" / {dtPnfdxdt=- ! {Vpn)\T)dx+ ! (Vp„)'(0)cix+ / / g^dtPndxdt. 

Jo J oj J (jj J (jj Jo J uj 

We estimate the RHS of the previous equation 

n{dtpnfdxdt < I I gndtpndxdt < \ j [ [dtpnfdxdt^-]- j j g^dxdt, 
Jo Joj ^ Jo Joj ^ Jo Joj 



and deduce j^{dtPnf'dxdt < J^^ J^g^dxdt. Then dtp G L^(0, T, L^(a;)) 
andatp-Ap = ^inL2(0,T,L2(a;)). □ 
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3.2 Proof of Theorem [23] 



We consider again the operator AC = -AC, D{A) = {C|C G H^{0),A( e 
L'^{0)} with its associated Dirichlet semigroup S{t). Thanks to Lemma 3.2, 
instead of considering directly the equation (12. 4p . in some of the following 
steps, we will consider the following equation 

u{t) = S{t){0) + [ f{u + v){s)ds= [ f{u + v){s)ds. (3.7) 
Jo Jo 

Step 1: We use a fixed point argument to prove that (13. 7p has a solution. 

We work with the case v G L°°(0, T, L'^p{0)). The case v E L°°(0, T, L°°(C)) 
can be treated with exactly the same manner. 

We consider the Banach space Yt = {u e L^^{{0,T), L'^p{0)),\\u\\yj, < 
oo} and ||^i||yj, = supo<t<Tt"||M(t)||2p- Let B be the closed ball in Yt with 
center and radius R2, we will use the Banach Fixed Point Theorem for the 
mapping 

$ : B^B (3.8) 

$(M)(t) = f S{t- s)f{u + v){s)ds. 
Jo 

Let ui, U2 be two functions in B, we now prove that $ is a contraction 

r||$(Mi)(t) - $(M2)(t)||2p < t° twSit - S)if{u, +V)~ fiu2 + V))\\2pds. 

Jo 

Using the — estimate (see Proposition 48.4 page 441 [21]), we obtain 

\\S{t-s){f{u, + v){s)~f{u2 + v){sm2,<{Mt-s)r^\\f{u,+v)-f{u2 + v)\\2. 

This leads to 

r t\\S{t - S){f{m + V)- f{u2 + V))\\2pds 

Jo 

< r(4^)-^ / it-s)-'-^\\f'iv + C)\\^\\ui-U2\\2pds 

Jo 

3(p-l) /■* 3(p-l) 

< t"(47r) 4p / (t-s) 4p ||(7^|^i^ + i + (7^|^i2 + \\^\\ui- U2\\2pds 

Jo "-^ 
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We have ||iti + ■i'||2p ^ (||i*i||2p + Iblbp)^""^, for i — Moreover, if p > 2, 



we have that (l^ilbp + Iblbp)^"^ < 2^-2(11^^11^^'+ \\v\\lp') and if p < 2, 

\2p ^ \ \^\\2p 



we have that (||Mj||2p+ |b||2p)^ ''" < ll'^illL^ + Il'^ll2p^) fo^ i — 1,2. These 



inequahties lead to 

\\S{t - S){f{u, + V){S) - f{u2 + V){sm2pds 



3(p-l) 

< r(47r)-^^C/max{2^'-^l} X 

t 



X I + 11^.2111,^ + 2||^;||^,-')|K-«2||2W. 



< r(47r)"T^C/max{2^'-2,l} / (t - 

Jo 



xis' 



3{p-l) 

< r(47r) 3^C/max{2f-2,l} X 



X / (t-60-^s-(^-^)"(2i?2 + 2r(f-i^'^||^||f;')||^/i-w2||2pcis 

< t°(47r)-T^C/max{2P-2,l}(2i?2 + 2T^-')"||^;||^-^p) x 

/"* 3(p-l) , 

X / {t-s)-^i^s~^-^^''\\ui-U2\\2pds 
Jo 

< t"(47r)-T^C/max{2f-^l}(2i?2 + 2T(f-^)°||^;||^^p) x 

/"* 3(p-l) 

X (t - s) 4p - ^2112^-5 P"(is 



< (47^)-^C/max{2^-^l}(2i^2 + 2^(f-l)1|^||L-^^)||«l-«2||y, x 

/"* 3(p-l) 

xt" / {t- sy^i^s^p^ds, 

Jo 

noticing that here we use the fact: ui,U2& B, i.e. [2^ |w2| |2p ^ 

< R2. 

We consider the last term on the RHS of the previous inequality 

* 3(P-1) 
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ft „, . „. 3(p-l) 



Jo 

3(p-l) 3(p-l) 

A simple computation gives 

/ 2 3(p-l) / 2 3(p-l) 
/ (1 - < / 4i^-P"c/l/ = 

JO -/o 



2 4p 



and 



/■■'■ 3(p-l) /■■'■ 3(p-l) 

/ (1 - uy^u-p'^du < / (1 - i^)-^-^"c/i/ 



3(P-1) 

As a consequence to these estimates 

r - + - /(«2 + V))\\2pds < 

Jo 

< (47^)-'^C^max{2^-^l}(2i^2 + 2T(^-^)"||^||^-^,) 



. 3(p-l) , 

= (47^)-'^Qmax{2^-^ l}(2i?2 + 2T(p''^^\\v\\'-^,^) x 

WeputCi = (47r)-^C/max{2^-^l}(2i?2+2T(^-^)ll^||L",2p)rS^^'^' 
thenCi < (47^)-^C^max{2^-^l}(2i?2+2Ti^-'^1|^;|e,)-^^^ < 



3(p-l) -, 3(p-l) p+3 

(47r)-— C/max{2^^-M}(2i?2 + 2||t;||^-^,)-^^^T/'' < i. This im- 
phes that $ is a contraction in the Banach space Yt- 
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Choosing U2 to be in this estimate, we obtain also that ||$(mi) — 
^(0)||yj, < Ci||Mi||yy < |i?2- Moreover, the estimate 

r f\\S{t-s){f{vm2pds < [\47r{t-s))-'^\\f{v)\\,ds 
Jo Jo 

< l\t-sr'V\\f{v)\Uds<\R, (3.9) 

implies that < i?2- Which means that $ is a contraction from a 

complete metric space to itself and it admits a unique fixed-point u in this 
set according to the Banach Theorem. 

Step 2: We prove that u G x (0,T)). 

Step 2.1: We consider the case v E L°^(0,T, L'^p{0)). 

Firstly, we prove that u e L°°{0,T,L^P{O)). From u{t,x) = S{t - 
s)f{u + v){s, x)ds, we have the following inequality 

\Ht,x)\\2p < f\\S{t-s)fiu + v)\\2pds< [\A7r{t-s))-'-^\\f{u + v)\\2ds 
Jo Jo 

< f\An{t-s)r-^{\\f{u + v)-f{v)\\2+mv)\\2)ds. 
Jo 

The computations 

\\f{v)\\2 < ||/(0)|b+||/VH|2< 11/(0)112+ ||/'(/i)||^|H|2p 

< ll/(0)||2 + ||C;|^;ri|U ||t;||2p=||/(0)||2 + C;||t;||^, 

p— 1 ^ 

implies that f{v) belongs to L°°(0, T, L^(C)). Combining this with the esti- 
mate of I |'u(t, |2p, we obtain 

\\u{x,t)\\2p < 

< fiinit - S)r-^{\\f{v + 011^ I \u\\2p + \\f{v) \U2)ds 

Jo "-^ 

< / (47r(t-s))-'^(||C^|^ + Cr^|U |M|2p+||/(t;)||oo,2)(is 

< j (47r(t-s))-^(||C^max{l,2f-2}(|^;ri + |Mri)||_^||M||2p 
+||/(^;)||oo,2)c^. 
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< j (4^(t-s))"^[C/max{l,2P-2 

+ \\f{v)\\oo,2]ds 



llp-1 I II I 



l2p^)ll^l|2p 



< 



Jo 



X 



X 



3(p-l)ii 
{4:71 {t - S)) 4p ds 

[\cfme.x{l,r-'}{\\v\f^;,^+\\u\\l^^^^^ 
Jo 



4p 



notice that h, I2 > 0, + = 1 and /i < 3(^;z:ij- 
The inequahty 



(t - s) *f (is = 



^ 3(p-l)ii * 
4p 



-, _ 3(p-i);i 

4p 



, 3(p-l)ii 



_ 3(p-i);i - _ 3{p-i)h 

4p ^ 4p 



, 3(p-l)ii 

T 5^ 



leads to 



4p 



3(p-l) 

(47r) 4p X 



(3.10) 



X 



/ [C^max{l,2^-^}(||^||L-i+||«||^;^)||«||2,+ ||/(t;)|U,2]'^rf^ 
Jo 



Put ^(t) = I I |2p; the inequahty fl3.10p becomes 
Uit)< 



I 3(p-l)ii X ^ 



T 4P \ ' 3(p-l)i2 

, 3(p-i)/. (4^)" X 

4p ' 



X 



r[C;max{l,2^-n(||^irJ,+ [/(s)'^)t/(s)^ + ||/(^;)|U,2r^ds. 
Jo 



!2. 

/ i_3(p-l)il \ i-^ 3(„_i)(„ 

We denote = f ^^^^J (47r)-"^ /J[C^ max{l, 2^-2} ( 

+ Uis)"^) U{s)^ +\\f{v)\\oo,2r'ds, then V'jt) = (^ ^"aSf ) " (47r)-^ 
[C; max{l,2*'-2| + j^^^)^) jj^^f^ + ||/(t;)||oo,2]'^ Since [/(t) < 



ip-i 

1 00, 2p 
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oo,2p 



V{t), then V'{t) < (47r)-^^^ [C^ max{l, 2^-^} ( 

+\/(.)'^) + ||/(t;)||oo,2]'^ 

Basing on the notations in (12. Sp . we put G{r) = 

3(p-l)i2 ,^ 

(47r) 4p 

/ ^ 3(p-l)ii \ ^ J 

[C^/max{l,2P-2}(||t;|r^J^ + C-)C^ + ||/(tOI|oo,2]'^ 

then G(r) is an increasing function on [0, +cx3) and G'{V{t)) < 1. Hence 
G{V{t)) < t + G{V{0)) = t + G{0) =t <T, and U{t) < V{t) < G-\T,). 

Finally, \\u{.,t)\\2p < G"^(r,)^. We have proved that u E L'^{0,T, L'^p{0)) 
by proving a New Gronwall Inequality on U{t). 

Secondly, we prove that u G L°°(0, T, L°°(C)). Using Proposition 48.4 
page 441, [21], we get the following estimate 



uit)\\^ < 



I oo 

< f\\S{t - s)f{u + v){s)\Uds < fiAnit - s)r^\f{u + v){s)\\2ds 
Jo Jo 







< I {A7r{t-s)rH\\f{u + v){s)-f{v)\\, + \\f{vmds 
{4n{t- s))-^Gfme^x{l,2P-'}{\\u\\^,;' +\\v\\^^^^^^^ 



< 







< r(47r(t-.))-i[Cymax{l,2P-2}(G-^(T,)^ + ||^;|r^-yG-^(T.)^ 
Jo 



+ \\fiv)\\oo,2]ds 

< ^-t(4r.)i[C;max{l,2^-2}(G'-i(T,)'^ + \\v\Li^)G-\T^f^ 



-wmw 



oo,2 



which deduces u G x (0,T)). 

Step 2.2: For the case v G L°°(0, T, we use exactly the same 

argument, but with the function G(r) = 

(47r) 4p dC 



/ 1 3(p-l)ii \ 

I 1 3(p-l)ii 
V 4p 



C^max{l,2P-2}((Mm(0)^)p-i) + C'^)C^ +max|^l<M|/(C)l"^(C)^ 
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Then we have \\u{t)\\^ < 7r-| (4T,)3 [Cf max{l, 2^-2} (G"! (T,) V +(Mm 

(C>)i)P-i) G-H^*)"^ +max|^|<M |/(C)I and m G L°°(C x (0,r)). 

Step 3: We prove that u G C{[0,T], L'^{0)) and m is also a solution of 
(12.41) . The proof in this step works well for both cases v G L°°(0, T, L'^^{0)) 
and v G L~(0,T,L~(C)). 

For e positive, 

+ e) - u{.,t) 

t+e j-t 

S{t + e-s)f{u + v){s)ds- / S{t- s)f{u + v){s)ds 

Jo 

t pt+e 

[S{t + t-s)-S{t- s)]f{u + v){s)ds + S{t + t- s)f{u + v){s)ds, 
which implies the following inequality 

I t + e)-u{.M2 < f\\ [Sit + e-s)-S{t~ s)]f{u + v) {s) W^ds + 

Jo 

+ ^ ^\\S{t + e-s)fiu + v){s)\\2ds. 

We firstly estimate the second term on the RHS of the previous inequality. 
Due to Proposition 48.4 page 441, 



[ ^\\S{t + e-s)fiu + v)is)\\2ds< [ ^\\f{u + v)is)\\2ds. 
Jt Jt 

Since J^~^^ \ \f{u + v){s)\\2ds tends to as e tends to 0, J^*^'' \\S{t + e — s)f{u + 
v){s)\\2ds tends to as e tends to 0. 

Now, we estimate the first term. Noticing that 1 1 [S(t + e — s) — S(t — 
s)]f{u + v){s)\\2 tends to as e tends to for all s in [0,t], and 

ll[S{t + e-s)-S{t-s)]f{u + v){s)\\2 
< \\Sit + e-s)f{u + v){s)\\2 + \\S{t-s)f{u + v){s)\\2<2\\f{u + v){s)\\2, 

where f{u + v) can be proved to belong to L°°{0, T, L'^{0)) by using the same 
argument as before; due to the Lebesgue Dominated Convergence Theorem, 
we deduce that 

lim / \\[S{t + e-s)-S{t-s)]f{u + v){s)\\2ds = 0. 
^-^0 Jo 
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Therefore \im,^o\\u{.,t + e) -v{.,t)\\2 = 0, or u e C{[0,T], L^{0)). 
Since v G C{[Q,T], L^{0)), u + v e C{[0,T], L'^{0)). Which imphes 
f{u + v) G C([0,T],L2(C)); and by Lemma u is also a solution of (Q- 

4 Proof of The Existence Results for Semilin- 
ear Heat Equation in a Cylindrical Bounded 
Domain - Theorem 12.21 



First, we consider the following equation 

dtv- Av = inQx (0,T,), 

V = g on (9fi, (4.1) 

v{., 0) = Mo(-) on fl. 

We recall the following definition: 

Definition 4.1. (Definition c page 69, ^) Let u be a domain in M^, we say 
that doj X (0, T] has the outside strong sphere property if for every Q = {xq, to) 
in duj X (0,T], there exists a hall K with center {x,i) such that K (1 {uj x 
[0, T]) = {Q} and \x — x\ > fi{Q) > 0, for every {x, t) inux [0, T], \t — to\ < e 
where e is a constant independent of Q. 

We will prove that, in fact dVt x (0, T] has the outside strong sphere prop- 
erty. 

Let Q = (xo,to) be on dfl x (0,T]. Since D is smooth, there exists x and a 
constant r > such that K = B^N{x,r) 0, = {xq} and \ \x — Xo||Kiv = r. 
We choose i = to and let {x',t') be in Bf^N+i{{x,t),r) n (0 x [0,T]), then 
\\x' — x||ir3 + \t' — t\ < r. Since \\x' — xll^iv < r; x' G Bf^N{x,r) fl and 
x' = Xq. We infer that t' = to and B^N+i{{x,t),r) n (fi x [0,r]) = {(xo,to)}. 
For {x, t) G r2 X [0, r], I |x — x\\^n > r > 0. Consequently, dQ x (0, T] has the 
outside strong sphere property. 

By Theorem 8 page 69, we can conclude that dQ x (0,T] has local bar- 
riers. Due to Theorem 5 page 123, [7], there exists a unique solution v to 
the problem (gl]), and v G C^''^{{0,T] x n) n C{n x [0,T]). According to 
Theorem 10 page 72, ^, v e x (0,T)). 
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Due to the Week Maximum Principle page 368, [H]: maXf^x[o,T] ^ = max:g(Qx(o,T)) 
V, and minf^x[o,T]^ = ming(nx(o,T)) ^- Which gives v < M and v > —M on 
n X [0,T] since M > max{| |iio| U, Hs'lloo}- Then \ \v\\L-^(^nx{o,T)) < M. 

Subtracting (14.11) and (12. 7p . and put w = u — v, we have the following 
equations 

dtw - Aw = f{w + v) infix(0,T), 

«;(.,.) = onaO, (4.2) 

u;(.,0) = onH. 

According to Theorem 12.11 (14.21) has a solution w satisfying w G L°°(f2 x 
(0,r))nC([0,T],L2(fi))nL2(0,r,/7oH^)) and e L2(0,T,L2(n)). More- 
over, ||w;||oo,oo<7r-| (4r,)3 [C/max{l,2P-2} (G'-i(r,)'i^ +((Mm(fi)^)P-i) 
G-\tS~^ +max|^l<M |/(C)|m(^])5], where G(r) = 

3(p-l)i2 

(47r) 4p 



3(p-l)ii \ I- 
rpl 3^ \ 1 



3(p-l)i 



4p 



C^max{l, 2P-2}((Mm(fi)^)P-i + C'^)C^ + max|^|<A/|/(C)|m(^])^ 



Hence (ETD has a solution u satisfying u e L°°(fix (0, T)) n ^([0, T], L2(fi)) 
nL2(0,r,i7i(n))andatnG L2(0,r,L2(fi)). Moreover, | U.oo < :=M+ 

TT-t (4r,)3 [C/max{ 1,2^-2} (G'-i(r,)'^ +((Mm(fi)i)P-i) G'-i(T,)'^ +max|^|<M 

Using the results in section VI. 8 [15] for the equation dtu — Au = g in 
Q X (0, T), where g = f{u + v), g & L°°{Q x (0, T)), we can conclude that u 
is continuous on f2 x (0, T). 
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5 Proof of the Well-Posedness and Conver- 
gence Properties of the Algorithms - The- 
orems 12.31 and 12.4 

5.1 The Well-Posedness of the Algorithm - Proof of 
Theorem [2731 

We consider the equation fl2.9p 



Since f{x) is positive for all x in M, then dt(f)M — ^4>ai > and dt{(pM — 
M) - A{(f)M - M) > . Since 0m = M on {dil x [0,ro]) U (fi n {0}), then 
<pM > M on Q X [0,Tq\ according to the Maximum Principle. 

We will prove that the algorithm is well posed for T < T* by recursion. 

For k = 1, we consider the problem on the j-th domain 



According to Theorem 12. 2[ the equation (15. 2p has a solution satisfy- 



ing u] G L~(fi,- X (o,T)) n c{[o,T],L\nj)) n L\o,T,H^{nj)), dtu] e 

L^(0, T, L^(fij))and is continuous on VLj x (0,T). The fact that e 
Ci'2((0,T) X fi)) then follows. 

Since dt{(l)M - - A(0a/ - li]) = /(0a/) - /(^^]); (?t(0Af - li]) - A(0m - 
mJ) - (0M - «J)c(a;,t) = 0, where 




(5.1) 



' d,u) - An] = /(«]) 

< m}(-,6j,-) = M°+i(-,f'i,-) 

^}(.,.)=^7(-,-) 
. w}(.,0)=wo(.) 



in fij X (0,T), 
in D X (0,r), 
in D X (0,r), 



on dD X [aj,foj] x [0,T], 
on f2j. 



(5.2) 
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We notice that c G C(^- x (0,r)) and u] < M < (pM on (5^- x [0,T]) U 
{flj X {0}). Then it follows by the Maximum Principle that (pM > ^} on 
X (0,T). 

Since ^((u) + M) - A(m] + M) - (/(^z]) - /(-M)) = /(-M) > and 
u] > —M on {d^j X [0, T]) U (fij x {0}); the Maximum Principle then implies 
that u] > -M on x (0,r). 

Consequently, (;/>m > u] > -M on x (0,r), Wj = Tj. 

Suppose that the algorithm is well posed up to step k = m, and moreover 
(pM > > —M on fij X (0, T) Vj = 1, /, z = 1, m, using the same argument 
as above we can see that the algorithm is well posed for k = m + 1 and 
0Af > > -M on X (0,T), Vj = TJ. 

Therefore, the algorithm is well posed for all A; G N and we have also 

> u'^ > -M on X (0,T), Vj = 



5.2 The Convergence of the Algorithm - Proof of The- 
orem 12.41 

We devide the proof into several steps. 

Step 1: The Exponential Decay Error Estimates. 

On the j-th domain, at the k-th iteration, we consider the equation 



r dte'; - ^e] = f\u]) - fiu) in X (0, T), 



I e,'(-,0)=0 



Clj , 



in D X (0,T), 
in D X (0,r). 
on dD X [aj,6j] x [0,T], 



(5.3) 



on 



Since the algorithm is well posed, according to Definition 12. ![ the set 
{u^\j = 1, /, A; G N} is bounded by a constant Ci. Which implies that {e^|j = 



1, /, G N} is bounded by a constant C2, where C2 = Ci + 1 1^| lc( nx(o r) )- 

We consider the Controlling Function $(a;, t) = P(e^) exp(/3(z — aj) — jt), 
and C{^) = dt^ — A$ + 2/35^$, where /3 and 7 are constants to be chosen 
later in the next steps. We now prove that if we choose (/?, 7) such that 
7 — is large enough, is negative. 
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A simple computation gives 

£($) = {dte'] - AeJ)P'(e5) exp{/3{z - a,) - -ft) + 

+ - 7)P(e,') expi^iz - a,) - ^t) - {We'^fP'\e]) exp(/3(;. - a,) - ^t). 
By (I5.3p . the previous computation leads to 

£($) = - f{u))P\e'^) exp(/3(z - a,) - 7^) + 

+ - 7)P(e5) exp(/3(z - a,-) - 7^) - {ye]fP'\e)) exp(/3(;. - a,) - 7^), 
which implies 

cm < [(/(«)) - f{u))P\e)) + - ^)P{e])] exp(/3(z - a,-) - 7^). 

We consider the term (/(-Uj) - f{u))P'{ej) + - 7)P(e^). By the Mean 
Value Theorem, f{uj{x,t)) — f{u{x,t)) = e'j{x,t)f'{(^{x,t)), where C(3^?^) 
is a number lying between u^{x,t) and u{x,t). We observe that < 
< Ci + 1 1^1 lc( nx(o,T) ) = ^'a, which implies \f'{C{x,t))\ is 
bounded on x (0, T). 
We fix a pair {x, t) in x (0, T): 

• If P(eJ(x,t)) = 0, then = and P'(e^^) = 0. That means (/(uj) - 
/(«))P'(et) + - 7)^(e?) = 0. 

• If P(e^^(x,t)) ^ 0, then 

(/(^^S - f{u))P'{e>;) + - 7)P(e^^) _ e';{x,t)P'{e';] 



P{e]) P{e]) 
Since is bounded by C2, ^'^^fj^k!^'^^^ is bounded by K{C2); which 



/(C(x,t)) + (/?2- 



means pr k-. ^ /'(C(^?^)) is bounded on i7 x (0,T) as |/'(C(x,t))| is 

bounded. Since P(e^) is positive, then if we choose 7) such that 
— + 7 is large enough, we can have that {f{Uj) — f{u))P'{ej) + — 
7)P(e^) is negative. 

Consequently, if 7 — /?^ is large enough, {f{uj±) — f{u))P'{e^) + {l3'^ — ^)P{e^j) 
is negative or is negative. 

When £($) < 0, according to the Maximum Principle, the function 
can only attain its maximum values on the boundary dVt x [0,T] 
or X {0}. We have that $ > 0, and moreover $ = on x {0} and on 
dD X [aj.hj] X [0,T]. Thus: 
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• If the maximum value(s) of $ can be achived on flj x {0} and on 

dD X [aj,bj] X [0,r], then $ = on Q x (0,r); which means that 
e'^ = on n X (0,T). 

• If 7^ 0, then the maximum value(s) of $ can be achived only on 
D X {aj} X [0, T] or on D X {bj} x [0, T]. 

For X in M''^, denote that x — {X, z) where X e R-'^"^ and z e M, we now 
have the following exponential decay estimates for the errors e^. 
Case 1: j = I. 

We have that ej^ = on D x {ai} x [0,T]. Since 7^ 0, the maximum 
valuc(s) of $ can only be achived on D x {61} x [0,T] and for {X,z,t) e 

fii X (o,r) 

P(et)(X,z,t)exp(/3(z-ai)-7t) < 
< max {P{e'l){X,bi,t) ex.-p{P{bi - ai) - -ft)} 

Dx[0,T] 

= max {P{et'){X, bi, t) exp(^S'i - 7^)}. (5.4) 

Dx[0,T] 



Case ^: j = I. 

We have that = on D x {bj} x [0,7]. Since e'} ^ 0, the maximum 
value(s) of $ can only be achived on D x {aj} x [0, T] and for {X, z, t) e 
0/ X (0, T) 

P(eJ)(X, 2;, t) exp(^(z - aj) - 7*) < 
< max {P(eJ)(X,a/,i) exp(-7f)} 

Dx[0,T] 

= max {P(eti)(X,a7,t)exp(-7t)}. (5.5) 

Dx[0,T] 



Case 3: 1 < j < I. 

The maximum value(s) of $ can be achived on both D x {aj} x [0, T] and 
D X X [0, T] and for (X, t) G 1^^- x (0,T) 

P(eJ)(X,2;,t)exp(/3(z-a,)-7t)< 
< max! max \P(e'^)(X, 6,-, t) exp(/3(6,- - a,) - jt)}, 

Dx[0,T] ■' 

max {P(e^')(X,a,-,t)exp(-7t)}} 

Dx[0,T] 
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= max{ max {P{e^_^_l){X, bj,t) exp{(3{bj — aj) — 7^)}, 

max {P(e^lJ)(X,a,-,t)exp(-7t)}}. (5.6) 

Dx[0,T] 

Step 2: Proof of convergence. 

Step 2.1: Estimate of the right boudaries of the sub-domains. 
Consider the /-th domain, at the /c-th step, (15. 5p infers 

P{e){X, z, t)) exp(/3(z - a/) - -ft) < max {P(eJ(X, a/, t)) exp(-7t)}. 

Dx[0,T] 

Replace z by we get 

P(eJ(X,6,_i,t))exp(/3(fe,,i - a,) -7t) < max {P(et(X, a,, t)) exp(-7t)}. 

Dx[0,T] 

Since ef(X,6,_i,t) = e^+|(X, t), 

P(e^+}(X,6,_i,t))exp(/3(6,_i-a,)-7t) < max {P(et(X, a,, t)) exp(-7t)}. 

Dx[0,T] 

Let (3 in this case be /?i to be a/I' 5 then when we choose 7 large, 7 — is 
large. The inequality becomes 

P(e,^+|(X,b,_i,t))exp(-7t) <exp(-/3i5,_i) max {P(et(X, a,, t)) exp(-7t)}. 

Dx[0,T] 

We deduce 

P{e'}tl{X,bj^ut))exp{-^t) < exp(-/3i5,_i)^fc. (5.7) 

Moreover, on the (/ — l)-th domain, at the {k + l)-th step, (15.61) leads to 

P{e'+l){X, z, t) exp(/?(^ - a,_i) - ^t) < 

max{ max {P(e^^l)(X, t) exp(/3(6/_i — a/_i) — 7^)}, 

Dx[0,T] 

max {P(e^+|)(X, a,_i, t) exp(-7t)}}. 

Dx[0,T] 

Since eJ+|(X,6,_2,t) = e'}tl{X,bj^2,t), 

P{e]tl){X, bj_2, t) exp(/3(fe,_2 - a/-i) - 7^) < 
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< max{ max {P(ej+|)(X, fo/.i, t) exp(/?(6/_i - ai^i) -7^)}, 
Dx[o,r] 

max {P{e]t\){X, ai.,,t) exp(-7t)}} 

Dx[0,T] 

= max{ max {P(e^+|)(X, fo/.i, t) exp(/?L7_i -7^)}, 

Dx[0,T] 

max {Pie'jtDiX, a,„i, t) exp(-7t)}}. 

Dx[0,T] 

Hence 

P{e'}tl){X,bj_2,t)exp{pSj.2-lt)<m8.x{P{e'}+l){X,bi_^,t) x 
X exp(/3L,_i - ^t),P{e'}+l){X, ai^i,t) exp(-7t)}. 

Since 

P(e^^J)(X,a,_i,t)exp(-7t) < E^+i, 

(15. 7p implies 

F(e^+2')(^>^/-2,i)exp(/357_2 - 7^) < max{^fc exp(/3L,_i - ^^+1}. 
Thus 

P(e^^2')(X,6,_2,t)exp(-7t) < max{Sfcexp(/3(L,_i-5,_2)-/3i>S/-i),^fc+iexp(-/35,_2)}. 

We choose (3 = (^2 = f^ij^ such that /32(— + S'/_2) + = /32Si_2] 

then 

F(e^^2')(^,^/-2,t)exp(-7t) < max{£;fc, exp(-/?25/-2). (5.8) 

Using the same techniques as the ones that we use to achive (15. 7p and 
(15.81) . we can prove that 

P{e^jtj){X, bi_j,t) exp(-7t) < max{£^fc, . . . , Ek+j-i} exp(-/3jS'/_j), (5.9) 

where /?, = /3ig^...|^,j = {2,. ..,/-!}. 

Step 2.2: Estimate of the left boundaries of the sub-domains 
Consider the 1-th domain, at the fc-th step, (15. 4p infers 

P(e'l(X, z, t)) exp(/?2 - 7t) < max {P(e\{X, 61, t)) exp(^6i - 7^)}. 

Dx[0,T] 
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Replace z by a2, we get 

P(et(X,a2,t))exp(-7t) < max {F(et(X, 61, t)) exp(-7t)} exp(/3(-a2+fei)). 

Dx[0,T] 

Since e^(X,a2,^) = e^+^(X, aa, t), 

P(e^+^(X,a2,t))exp(-7t) < max {P(e?(X, 61, t)) exp(-7t)} exp(^5i). 

Dx[0,T] 

Let = —/3[ = —\f\-, then when we choose 7 large, 7 — is large. The 
inequality becomes 

P(e2'+i(X,a2,t))exp(-7t) <exp(-/3;5i) max {P(et(X, 61, t)) exp(-7t)}. 

Z)x[0,T] 

We deduce 

F(e^-+i(X,a2,t))exp(-7t) < exp(-/3;>Si)^fc. (5.10) 
Moreover, on the 2-th domain, at the ik + l)-th step, (15. 6p leads to 

P(e^+^)(X,z,t)exp(/3(;2-a2) -7^) < 

max{ max {F(e2"*'^)(X, 62, t) exp(/3(62 -02) -7O}) 

Dx[0,T] 

max {P(e^+^)(X,a2,t)exp(-7t)}}. 

Since e^+^(X, 03, t) = eg+^(X, 03, t), 

F(e^+2)(X, a3, t) exp(/3(o3 - as) - 7^) < 

< max{ max {P[eX^''){X, 62, t) exp(/3(62 - 02) - 7^)}, 
Dx[o,r] 

max {P(e^+^)(X,a2,t)exp(-7t)}} 

Dx[0,T] 

= max{ max {P(e2+^)(X, 62, ^) exp(/3L2 - 7^)}, 

Dx[0,T] 

max {P(e^+^)(X,a2,t)exp(-7t)}}. 

Dx[0,T] 

Hence 

F(e^+2)(X,a3,t)exp(/3(L2->S2) -7O < max{P(e^+^)(X, 62, t) x 
X exp(/3L2 - 7t), /'(e^+')(X, as, t) exp(-7t)}. 
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Since 

P(e^+i)(X,62,t)exp(-7t)<i?,+i, 

(15. 7p implies 



P(e3^+2)(X, aa, t) exp(/3(L2-52)-7i) < max{^fc exp(-/3;>Si), ^fe+i exp(/3L2)}. 
Thus 

P{et^^){X, as, t) exp(-7t) < max{^fc exp(-/3(L2-52)-/3i5i), ^fc+i exp(/352)}. 

We choose /3 = = -^^g such that -/3(L2 - ^2) - = PS2. Then, 
we have that 

F(eg+2)(Xa3,t)exp(-7t) < max{^fc, ^^+1} exp(-/3^52). (5.11) 

Using the same techniques as the ones that we use to achive (15.101) and 
(15. lip , we can prove that 

F(eJ+^"^)(X,aj,t)exp(-7t) < max{£'fc, . . . .,Ek+j-2} e^'?{-P'j^iSj-i), 

(5.12) 

where = g . . . j = {2, . . . , J - 1}. 
Step 2.3: Convergence resuh 
From the fact that e = \/? f^"'f^~^ , and 

Ek = max {£^A:+j}, 

and (15321) give us 

Efc+i <Efcexp(-e),VA;eN, 

or 

< ^oexp(-ne), Vn G N 

Hence i^^ tends to as /c tends to infinity. Which gives 
^lim mg||P(ef)||c(^^-^) = 0, 

and we get the linear convergence of the algorithm. 
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